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Abstract. We solve the stationary Navier-Stokes equations for 
non-Newtonian incompressible fluids with shear dependent vis- 
costy in domains with unbounded outlets, in the case of shear 
thickening viscosity, i.e. the viscosity fi is given by the power law 
/z = |D(v)|P~^, where |-D(v)| is the shear rate and p > 2. The flux 
assumes arbitrary given values and the Dirichlet integral of the ve- 
locity field grows at most linearly in the outlets of the domain. Un- 
der some smallness conditions on the "energy dispersion" we also 
show that the solution of this problem is unique. Our results are 
an extension of those obtained in [15] for Newtonian fluids {p = 2). 



1. Introduction 

The Navier-Stokes system for stationary incompressible flows in 
a domain with unbounded straight outlets, with the velocity field con- 
verging to parallel flows (Poiseuille flow) in the ends of the outlets, was 
solved first by C. Amick |2] in the 1970s. This problem is known as 
Leray problem, cf. [21 P- 476]. Amick's solution assumes the fluxes of 
the fluid in the outlets to be sufficiently small, which turns out to be 
a sufficient condition to deal with the convective (nonlinear) term in 
Navier-Stokes equations. It is an open problem to solve Leray problem 
for arbitrary fluxes. Alternately, Ladyzhenskaya and Solonnikov [T5] 
considered the stationary Navier-Stokes equations not demanding the 
fluid to be parallel in the ends of the outlets, but instead having arbi- 
trary fluxes. In this case, the outlets do not need to be straight and 
they solved this new problem for domains having arbitrary uniformly 
bounded cross sections and with the fluid having arbitrary fluxes. Be- 
sides, their solution has the property that the Dirichlet 's integral of 
the velocity field of the fluid grows at most linearly with the direction 
of each outlet, and they also proved that this solution is unique under 
some additional smallness condition. 
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In this paper we extend the Ladyzhenskaya-Solonnikov's theorem, 
i.e. "Theorem 3.1" in [15], for power- law shear thickening fluids, i.e. 
incompressible non-Newtonian fluids obeying the power law 

(1.1) § = |D(v)r^D(v), 

when p > 2. Here, § is the viscous stress tensor, v is the velocity field 
of the fluid and -D(v) is the symmetric part of velocity gradient Vv 

(i.e. A,(v) = Kll + g-) for V = (t;i,---,^;0, ^, J G {1, ■ ■ ■ , n}, 
n G N). For p = 2, the fluid is Newtonian. If 1 < p < 2, the fluid is 
called shear thinning (or plastic and pseudo-plastic) and if p > 2, shear 
thickening (or dilatant). In engineering literature the power law f ll.ip 
is also known as Ostwald-De Waele law (see e.g. [B]). Corresponding 
to f ll.ip we have the following system of equations modelling the flow 
of an incompressible fluid in a stationary regime: 

r -div(|D(v)r2z)(v)) + (vV)v + VP = 
^ ' \ divv = 0, 

where V is the pression function of the fluid (and v is the velocity field, 
as already indicated above). This model equations are also referred to 
as Smagorinsky model, due to |22j, or Ladyzhenskaya model, due to 
[T^ [T5| [13]. A related model where the viscosity is given by |vp~^, 
instead of |D(v)|p~^, is considered in |16j. For this case, it is shown in 
[T6t Remark 5.5 in Chap. 2, §5.2] the existence of a (weak) solution for 
system (11. 2p in a bounded domain with homogeneous Dirichlet bound- 
ary condition, for p > . There are many results concerning the 
solution of (11. 2p in bounded domains. For instance, in [9] the exis- 
tence of a solution for (11. 2p is obtained under the weaker condition 
thatp>;^. 

In unbounded domains there are not so many results. For parallel 
fluids we can identify v with a scalar function v and the system (II. 2p 
reduces to the p-Laplacian equation 

(1.3) - diY{\Vv\P-^Vv) = c 

for some constant c (related to the ^^pressure drop"). So, we can con- 
sider the Leray problem for (II. 2p . i.e. the solution of (II. 2p in a domain 
with straight outlets with the velocity field tending to the solution 
of (11. 3p in the ends of the outlets. This problem was solved by E. 
Marusic-Paloka [T7] under the condition that the fluxes are sufficiently 
small and p > 2, thus extending Amick's theorem [2] for power fluids 
with p > 2. As far as we know, the Leray problem for (II. 2p when p < 2 
(with small fluxes) is an open problem. 

In this paper, as we mentioned above, we extend Ladyzhenskaya- 
Solonnikov's theorem [T5j Theorem 3.1] for (II. 2p when p > 2. More 
precisely, we obtain the existence of a solution v to the system (II. 2p for 
n = 2, 3, and p > 2, in a domain Q with unbounded outlets, specified 
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in the next Section, for any given fluxes in the outlets and homoge- 
neous Dirichlet boundary condition vldQ = 0. The "Dirichlet inte- 
grals" J I Vv|P of our solution grows at most linearly with the direction 
of the outlets (see (12.11) ^ in Section H]). Besides, we observe that these 
integrals over portions of the outlets with a fixed 'length' are bounded 
by a constant that tends to zero with the flux (see Proposition 14.21 
and Remark H]). Under this condition and some aditional one, we have 
uniqueness of solution (see Theorem 14. 4p . All these facts were obtained 
in [15] for the case p = 2, but the power-law model ( fll.2p with p ^ 2) 
was not treated in [TJ]. In the next two paragraphs we look at some 
facts relating to the case p 2. 

First, to deal with the nonlinear term div(|D(v)|''~^D(v)) one can 
use the monotone method of Browder-Minty. Secondly, we extend the 
technique employed in [15] to obtain the existence of a solution, which, 
in particular, consists in first solving the problem in a bounded trun- 
cated domain and then taking the limit when the parameter of the 
truncation tends to infinity, to obtain a solution in the whole domain. 
To take this limit we need first uniform estimates with respect to the 
truncation parameter for the solution in the truncated domain, and 
this is obtained by integrating by parts the equation times the solution 
in some fixed bounded domain. Then we need the regularity of the 
solution in bounded domains, more precisely, that the solutions have 
velocity field at least in the Sobolev space W^''' and pressure in W^'\ 
for some positive number I, due to the boundary terms that comes from 
the integration by parts. However this regularity is not expected for 
the weak solutions of (11. 2p . if p 7^ 2. To overcome this difficult, when 
dealing with (11. 2p in a truncated bounded domain we modify it to 

r -div{(i + |D(v)r2)D(v)} + (v.v)v + vp = o 

\ divv = 0, 

where T > is the truncation parameter. See Proposition 14.11 in Sec- 
tion H 

As in [23] and [15] , and in several subsequent papers, here the velocity 
field V is sought in the form v = u-|-a, where u is the new unknown with 
zero flux and a is a constructed vector field carrying the given fluxes in 
the outlets (i.e. if the given flux in an outlet with cross section S is a 
then a-n = a and u-n = 0, where n is the unit normal vector to S 
pointing toward infinity). This vector field a depends on the geometry 
of the domain and, in the aforementioned papers, its construction is 
very tricky and makes use of the Hopf cutoff function (see [231 fT5]). 
In the case of power-law fluids (II. 2p with j9 > 2 we found out that 
the construction of a can be quite simplified. Indeed, a key point in 
the construction, in any case, is to obtain a vector field a that controls 
the quadratic nonlinear term (uVu)a, which appears after substituting 
V = u + a in (II. 2p and multiplying it by u. That is, to obtain a priori 
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estimates, one multiplies the first equation in fll.2p by u and try to 
bound all the resulting terms by the 'leading' term |iI)(u)|P. In [TB] it 
is shown that for any positive number S there is a vector field a which, 
in particular, satisfies the estimate 

|una|'<c5' / |Vu|^ 
nt Jnt 

for some constant c indepedent of 6, u and fit, where Qt is any trun- 
caded portion of the domain with a length of order t. Looking at their 
construction and using Korn's inequality it is possible to show that 



1.5) / \uf\sLf <c6P't^P-^^/^P-^^ ( [ \D{u)A 

Jnt \Jnt J 



v' Iv 



where p' is the conjugate exponent of i.e. p' = p/{p — 1). When 
p = 2 this estimate reduces to | /^^(uVu)a| < c6 J^^ |Vup. With 
this estimate we can estimate the integral of (uVu)a in the truncated 
domain Qt, by using Holder inequality: 

1/2 / /• \ 1/2' 

2 1 / / l„|2'|„|2' 



1.6) 



(uVu)a| < / |Vu|2 / |u|2'|a 



nt 

< cS I |Vu|2. 

nt 



Thus we can control the nonlinear term (uVu)a by taking necessarily 
S sufficiently small. When p > 2, proceeding similarly and using also 
Korn's inequality, we obtain 



(1.7) I [ (uVu)a| < c(5t(^'-2)/P ( [ \D{u)A 

J nt \J nt J 

Then, by Young inequality with e, we have 



2/p 



(uVu)a|<e/ \D{u)\P + C,t, 
nt J nt 

for some new constant C^. From this estimate, we can control the 
nonlinear term (uVu)a by taking e sufficiently small, and so we do 
not need to construct the vector field a satisfying the estimate (11. 5p 
for a sufficiently small S. See Section H] for the details. In fact, if a is 
only a (smooth) bounded divergence free vector field vanishing on dQ, 
then, by Poincare, Holder and Korn inequalities, and the fact that our 
domain has uniformly bounded cross sections and p/p' = p — 1 > 1 
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{p > 2), we have 



\u\P \a\P < c I \u\P < c / |Vu|P 

Ut J Ut 

p'/p 

< ct^-P'/P i I \Vu\P 



,1.8) . ^p'/p 

= ct^P-2)l{p-l) / IVulP 



< ct{p-2)/(p-i) f / \D{vi)\P '''''' 



which is (11.51] for 5 = 1. 

The plan of this paper is the following. Besides this introduction, in 
Section[2]we introduce the main notations and set precisely the problem 
we will solve, state a lemma about the existence of the vector field a, 
carrying the flux of the fluid, and state our main theorem (Theorem 
I2.2|) . In Section [3] we state some preliminaries results we need to prove 
our main results. In Section H] we prove our main theorem, make some 
remarks and prove a result about the uniqueness of our solution. 

2. Ladyzhenskaya-Solonnikov problem for power-law 

FLUIDS 

In this section we set notations and the problem we are concerned 
with and state a lemma and our main theorem. 

We denote by a domain in R", n = 2, 3, with a boundary, of 
the following type: 

2 
1=0 

where Qo is a bounded subset of M", while, in different cartesian coor- 
dinate system, 

Qi = {x = {xi,x') G M"; Xi < 0, x' G Si(xi)} 

and 

Q2 = {x = ixi,x') G M"; Xi > 0, x' G ^2{xi)}, 

with Si(xi) being C°° simply connected domains (open sets) in R"~^, 
and such that, for constants /i, /2, < /i < ^2 < C)0, they sastify 

sup diamSj(xi) < I2 

(-l)'xi>0 

and contain the cylinders 

CI ={xe M"; i-iyxi > e \x'\ < |} C ^ ; 
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For simplicity, we will denote by S any of the cross sections Sj = 
Ej(xi) or, more generaly, any cross section of Q, i.e., any bounded 
intersection of fl with a — -dimensional plane. Wc will denote by 
n, the ortonormal vector to E pointing from Qi toward Q2 i-e. in the 
above local coordinate systems, we have n = (1, 0) (where e M"~^) 
in both outlets Qi and fl2- With these notations, the flux through any 
cross section E of Q of an incompressible fluid in fl with velocity field 
V vanishing on dQ, is given by the 'surface' integral J^, v ■ n (notice 
that by the divergence theorem applied to the region bounded by dQ, 
El and E2, we have J^,^ v ■ n = v ■ n, for any cross sections Ei and 
E2 of ill and ^2, respectively). 

We remark that we take our domain Q with only two outlets fli, 
i — 1,2, just to simplify the presentation, i.e. we can take Q with any 
finite number of outlets with no significant change in the notations, 
results and proofs given in this paper. 

We shall use the further notations, where U is an arbitrary subdo- 
main of fl, s > t > and 1 < g < 00: 



Qi^t = {xe Qi^i-iyxi <t},i = 1,2 



i,t,s — ^ ^i,5 \ ^ H,t 
Qt.s = Qs ^ Qt 



= (/Jvr+|Vv|'') 



1/9 



|v|l,.,^=(/t;|Vv|«) 

(u,v) = (u,v)f^ = Jf/U- V 

V{U) ^{cpe C^{U; R"); V • = 0} 

P^'^(C/) = ©(CT)'''''" 



In these notations, the set f2j - a bounded cut of with a "length" 
of order t - will be taken usually for large t, so this notation will not 
cause confusion with the (unbounded) outlets Jlj, where i — 1,2. 

By W^''i{U) and W^'^iU) we stand for the usual Sobolev spaces, 
consisting of vector or scalar valued functions, and Wfl^^{U) is the set 
of functions in W^''^{V) for any bounded open set V G U. Often when 
it is clear from the context we will omit the domain of integration in 
the notations. 

The notation \E\ will stand for the Lebesgue measure of a Lcbcsgue 
measurable set E in the dimension which is clear in the context. Finally, 
the same symbol C, c, C. or c. will denote many different constants. 

In this paper, we are concerned with the following problem: given 
q; e IR, find a velocity field v and a pressure V such that 
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(2.1) 



' div{|D(v)|P-2/}(v)} = vVv + VP inn 

V-v = inQ 

V = on dVl 

/^v^n = a 
supj>ot < oo. 

Cf. Problem 2.1 in (for the case p = 2). Here, and throughout, we 
use the notation 

V ■ Vw = (v ■ V)w = ^ Vi- — w = (v ■ Vwi, ■ ■ ■ , V ■ VWn) 

OXi 



i=l 



for any velocity fields v = (t>i, ■ ■ ■ , t>„) and w = {wi, ■ ■ ■ , Wn) defined 
in Q such that the last expression on the right makes sense. 



To solve (12. ip . we seek a velocity field v in the form v = u + a, 
where u is a vector field with zero flux and a will carry the flux a, i.e. 

u ■ n = and a ■ n = a. More precisely, we shall take a to be a 
vector field having the properties given by the following lemma. 

Lemma 2.1. For any p > 2, there exists a smooth divergence free vec- 
tor field a, which is bounded and has bounded derivatives in Q, vanishes 
on dQ, and has flux one, i.e. Jj, a = 1 over any cross section E of fl. 
In particular, given a G M, the vector field a = aa is a vector field 
preserving all these properties but having flux a and else satisfying the 
following estimates: 

i) 4|aKl^K<c|ar>-2)/(^-i)|^|?;^,^^, Vt>0, \f^eV{n); 
io. t-i t I'^^l^ - Vt > 1, i = 1, 2; 

iii) /f^jVap < c\af{t + 1) , Vt > 1; 
where p' = p/ (p — l) and c is a constante depending only on a, p and Q. 

The proof of this lemma is given in Section HJ 
For V = u + a, problem (12. ip becomes 



(2.2) 



f dw{\D{u) + D{a)\P-\D{u) + D{a))} 

= u- Vu + u- Va+a- Vu + a- Va + VP inn 
V ■ u = inQ 
u = on d^l 
k u ■ n = 

Formally, multiplying (I2.2p i hj (p = {ipi, ■ ■ ■ , </?„) G V(Q) and noticing 
that 



(2.3) 



'^dxj . . 
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and W ■ V = ~ ■ y = 0) after integration by parts we get 

|D(u) + D{^)r\D{u) + D(a)) : D{^) = 
^ ■ ' - (u ■ Vu, v?) - (u ■ Va, c^) - (a ■ Vu, c^) - (a • Va, cp) , 

for all (p G "DlQ), where for n x n matrices A = (aij), B = (hij) we use 
the notation A : B = 'Y^^j=i O'ijbij- Thus, the following definition for a 
weak solution to the problem fl2.2p is in order. 

Definition 1. A vector field u is said to be a weak solution to the 
problem fl2.2p if it has the following properties: 

i) u G wtm; 

ii) u satisfies fl2.4p for every cp G T>{Q); 
iii) u satisfies (1221)2-(I221)5- 

Similarly, a vector field v is said to be a weak solution to the problem 
(EH) z/ V G anc? satisfies (EI]) 2 -([21]) 5 arw^ ([221) wt/i u + a 

replaced by v, z.e. 

/ |D(v)r2Z}(v):I)(<^) = -(vVv,<^) 
/or a// (p G ^'(^]). 

Remark 1. The use of divergence free test functions cp in fl2.4l) e/im- 
inates the pressure V, but it is a standard fact that it can be recovered 
due to 'De Rham's lemma' (cf. e.g. [10, Lemma IV.1.1]^. 

We end this Section stating our main theorem, which we prove in 
Section HI 

Theorem 2.2. Let p > 2. Then, for any a G R, problem (12. ip has a 
weak solution v, in the sense of DefinitionUl 

3. Preliminary results 

In this Section we give some preliminary facts we shall need to 
prove our mains results in Section |H We begin with Lemma [3. II below, 
which is due to Ladyzhenskaya and Solonnikov |15l Lemma 2.3]. Our 
statement below differs slightly from [15] and, for convenience of the 
reader, we present its proof, which essentially can be found in [15] and 
[201 EI]. 

Lemma 3.1. Let \E' : M -^■ R 6e a strictly increasing function, 6 be a 
number in the interval (0, 1) and to < T . 

i) If z and ip are difjerentiable functions in the interval [to,T] satis- 
fying the inequalities 

r z{t)<^{z'{t)) + {i-5Mt), 

for all t G [to, one? z{T) < ip{T), then 

z{t)<^{t), VtG [to,T]. 
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ii) Suppose \I/(0) = 0. If z : [to, oo) [0, oo) is a non identically zero 
and non decreasing differentiable function, and satisfies the inequality 
z(t) < "^{z'it)), for all t > to, then limt_^oo ^(^) = oo. Besides, if 
^(r) < cr™ for all t > Ti, for some constants m > l,c > and 
Ti > 0, then 

— ^ 

lim inf t '"-i z{t) > ; 
If, however, \&(r) < or, for t > ti, then 

lim inf e~^^''z{t) > . 

Proof, i) Suppose that (p{ti) < z{ti) for some ti G [to,T). Then, 
by the first inequahty, we have -2(^1) < 6'^"^ {z' (ti)) , and so, using 
the second inequahty, we have also 6'^"^ {ip' (ti)) < (p{ti) < z{ti) < 
5-^^(2;' (ti)), then, < Since is strictly increasing, 

it follows that z'{ti) > f'iti). Consequently, z{t) > ip(t) for all t on a 
neighborhood on the right of ti, and so, taking t2 to be the supremum 
of these points in (ti,T), we have ti < t2 < T and, by the previous 
reasoning, we have z'{t) > f'{t) for all t in {ti,t2), but this yields a 
contradiction, since z{t) — (p{t) > is strictly positive ai t = ti and 
must be zero at t = ^2- 

ii) Let ti > to such that z{ti) > and A = '^^^{z{ti)). Notice 
that A > 0, since ^E'(O) = and is strictly increasing. As z is a 
nondecreasing function, we have that z{t) > z{ti) for all t > ti. Then 
we claim that z{t) > z{ti) + X{t — ti) for all t > ti. Indeed, the 
inequalities z{t) > z{ti) and z{t) < ^ imply z'{t) > ^-^(^(t)) > 

\E'~^(z(ti)) = A. Thus, we have shown the first statement in part 2) of 
the Lemma. For the remainder, notice that, since limt_!.oo z{t) = 00, 
there exists a r such that z{t) > ti for all t > r, so from "^{t) < cr"^ 
and z{t) < 'i'{z'{t)) we have z{t) < c(z'(t))'" for all t > r, and the 
results then follow by direct integrating this inequality. □ 

In the next lemma we collect three very useful inequalities. The first 
can be found in many texts, as for instance in [7] and [5], Lemma 2.1, 
p. 526]. The third inequality contains Korn's inequality (see |18j*l ). 
The last one is a classical Poincare type inequality; see e.g. [TOl p. 56]. 
In these inequalities, ci,C2 are positive constants depending only on p 
and, for the last two, on the domain U. 

Lemma 3.2. 

i) 

(|a;|P"^a; - \y\P~'^y,x- y) > ci\x - y\'^ + \y\^''^) > C2\x - y\P 

for all X, y G M" and p > 2. 

*In [TB], Korn's inequality is stated for dimension three. The resuh in dimension 
two can be obtained from the one in dimension three by extending the domain 
f/ C R2 to [/ X (0, 1) and the vector field v : [/ ^ ^o (v, 0) : [/ x (0, 1) ^ M^^ 
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ii) 

Cl|v|i,p,(7 <|| D{v) \\p^u< C2\^\l,p,U, 

for all V G V^'P{U) such that v|r = 0. 
iii) 

II V \\g^u< ci {\v\i,q,u+ II V ||i,r) , 
for all V G W^''^{U). In ii) and iii), U is an arbitrary bounded domain 
o/M"", n = 2,3, with a smooth boundary, T is any Lebesgue measurable 
subset of dU with positive measure, and 1 < p < oo. 

Next, we state a corollary of Brouwer fixed point theorem. 

Lemma 3.3. Let F : R" — > R" be a continuous map such that for 
some p > 0, ■ ^ ^ for all ^ G R" with \^\ = p. Then, there is a 
^0 e M" with \^o\<p such that F(^o) = 0. 

For a proof, see |IDl Lemma VIII.3.1] or ^ p. 493]. 

The next lemma yields a solution v of the equation div v = / satis- 
fying a nice estimate. This result is an important step in the proof of 
our main theorem. 

Lemma 3.4. Let U be a locally Lipschtzian and bounded domain in 
M", n > 2, and 1 < g < oo. Then there is a constant c such that, for 
any f G L'^{U) satisfying f = 0, there is a vector field v G Wq''^{U) 
such that V ■ V = / and \\ v ||i,g,(7< c || / \\g^u. 

See [ini Teorema III.3.2]. 

The final result of this Section regards the regularised distance func- 
tion to the boundary of a domain (an open connected set) in M". 

Lemma 3.5. Let V be a domain in M" and d{x) = dist{x, dV), x & V. 
Then, there is a function p G C°°{V) such that for every x V and 
any derivative , (3 = ■ ■ ■ , /3„) G Z+, we have 

, . d{x) < p{x) and 

^ ' \d^p{x)\<kp{d{x)Y-\^\ 

where kp is a constant depending only on (3 and n. 
See Theorem VI. 2]. 

4. Proof of Theorem 12.21 and other results 

In this section we prove Lemma 12.11 and our main theorem - The- 
orem 12.21 Besides, we make some remarks, prove a Proposition on 
the 'uniform' distribution of energy dissipation (Proposition 14.21) and 
a Theorem regarding the uniqueness of solution of problem fl2.ip . 

We begin by proving Lemma [2?T1 As we observed in the Introduction, 
the proof of this lemma (the construction of a) is simpler in this paper 
(i.e. for the case p > 2) than for the classical one for newtonian fluids 
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{p = 2). For the construction in the case p = 2, see ^15jj p. 744] and 
references therein; see also jTHl Lemma XI. 7.1, p. 272] and [121 p. 46]. 

Proof of Lemma 12. IL Suppose we have a vector field a as in Lemma 
12.11 Then the statements with respect to a = aa follow, with c de- 
pending on p, sup|^.^|>Q supq |a| and supj^ |Va|. Indeed, for property 
Lemma [2. H i), see (11.81) . For property ii), we have 

jVa|P < (sup |Van(sup |S|)|ar 

and iii) follows from ii): 

LjVal^ = LolVap + E.=i,2 4JVar 

< |fio|(sup |Van|a|^ + ((sup |Va|P)(sup |S|)) \aft. 

To construct a vector field a with the properties in the statement of 
Lemma [2?n first we observe that it is enough to construct in each outlet 
fli a vector field a* satisfying these properties in fij. Indeed, if we have 
this, then we can obtain the desired vector field a defined in Q by using 
appropriate cutoff functions. We omit this part of the proof and refer 
to [TOl cap. VI] for a similar procedure in a domain with straight outlets 
and Poiseuille flows in place of the vector flelds a*, to be constructed 
below. 

We flrst construct a in the case n = 2. By what we observed above, it 
is enough to construct the vector fleld a in an arbitrary outlet Qi , which 
we shall denote by Q in this proof. Without loss of generality, we take 
Q = {x = {xi,X2) G ; /i(xi) < X2 < /2(a^i)} for smooth functions 
/i,/2 such that fi{xi) < | < /2(xi) and f2{xi) - /i(xi) < k, 

for all Xi G M. {h < I2 are positive numbers introduced in Section O) 
Then we set 

a = V\ = (9.,C,-5x,C), 
for C(2^i5 X2) = ip{x2/ p{x)), where p(x) is the regularised distance to dVt 
(see Lemma 13.51) and t/^ : R — )■ M is a smooth nondecreasing function 
such that ip{s) = if s < and 1, if s > 1. We notice that C is identi- 
cally zero in the 'lower band' {x G f2 ; /i(xi) < X2 < 0} and identically 
one in a neighborhood of the 'upper boundary' {x G dVL ; X2 = f2{xi)}. 
In particular, a is a divergence free bounded vector fleld vanishing on 
a neighborhood of dQ and 

/ a-n= / Cx2dx2 = C{xi, f2{xi)) - C{xi, fi{xi)) = 1 . 
Jt. Jt, 

Now, because C, is constant in a neighborhood of each of the two com- 
ponents of dVt, we have that any derivative of C, is zero in this neighbor- 
hood and, thus, bounded in Vt. Then a and its derivatives are bonded 
function in VL. 

In the case n = 3, we take C,{xi^x') = ip^x']/ p{x)), x' = {x2,xs) G 
M^, where p{x) is the regularised distance to dQ (see Lemma [3.51) . ip is 
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(-X3,X2). 



Notice that ( constant for x' close to zero and equal to one in a neigh- 
borhood of dQ (i.e. p(x) close to zero), and thus, C is a smooth function 
with bounded derivatives, vanishing in neighborhoods of x' = and dQ. 
Therefore, a is a smooth function vector with bounded derivatives. Be- 
side, it is divergence free, and, by Stokes theorem in the plane, we have 
Ie^-^ = Ids^da = 1. □ 

To solve problem (12. 2p . first we shall solve the truncated modified 
problem, T > 0: 



Then we will use Lemma 13.11 to obtain a weak solution of (12. 2 p by 
taking the limit, when T — t- oo, in the solution u-^ of (14. ip . extended 
by zero outside Qt- 

Proposition 4.1. Let p > 2 and T > 0. Then problem (14. ip has a 
solution (u^,P) m vI'^{Qt) x Lp{Qt) n W^^^Qt) x W^^^Qt), for any 
t e (0, T), where I = 2q/ {p + q — 2), being q = 2p + 2 if n = 3 and any 
number in [1, oo) if n = 2. 

Proof. The regularity part, i.e. (u^,^) G W^'\Qt) x W^'\Qt), for any 
t G (0, T), is a corollary of the proof of Theorem 1.2 in [1]. Notice that 
if (u^, V) is a weak solution with in Pq'^ then v = u-^ + a is a weak 
solution in W^'P{nT) of 



r div{(i + |D(v)|P-2)D(v)} + v Vv + VP = infi. 
(4.2) < V • V = in fir 
I V = a on dilx- 



The fact that we do not have here the homogeneous Dirichlet boundary 
condition v = here in the whole boundary OQt does not affect the 
method given in ^ because a = in (Sfi^) H (dQ) and the remaing 
part of dQx, i-e., (dflT) / (dfl) , is interior to Qt- 

Then we have only to show the existence of a weak solution for (14. ip . 
For simplicity, most of the time in this proof we shall write Qt = ^ 
and = u. Also we keep the notation (-, ■) with the integration 
over Q = in this proof, i.e. for (vector) functions v, w such that 
v-w e L^lff"), (v, w) = J^j, v-w. We will apply the Galerkin method 
and the monotonicity method of Browder-Minty (cf. ^ Remark, p. 
497]). The Browder-Minty method is used due to the nonlinear term 
in the left hand side of (I4.ip i . 



(4.1) 
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Let {f^'jj = 1, 2, ■ ■ ■ } C V{Q) be a denumerable and linearly inde- 
pendent set of functions whose linear hull is dense in We shall 
write for m = 1, 2, . . ., 

m 

(4.3) u™ = 5^c,„,yp^ 

i=i 

where (cim, ■ • • , Cmm) £ solves the algebraic system 
(4.4) 

+ |D(u'") + D{a)\P-\D{n^) + Z}(a)) : D{ip^) 
+ (u"^ ■ Vu™, ip^) + (u™ ■ Va, ip^) + (a ■ Vu'", ip^) + (a ■ Va, ip^) = 0, 

j = 1, . . . ,m. To see that system (14.41) has a solution (ci^, ■ ■ ■ , Cmm), 
let F = (Fi,--- ,F„) : M'" ^ be the map such that, for ^ = 
(6, • • • ,^m) e M™, Fj(0 is defined by the left hand side of (jOD with 
_ 'Y^^iij^^ ■ By Lemma [3.31 it is enough to show that there is a 
p > such that F(0 ■ ^ > for all |^| = p. Since {vT ■ VvT, vT) = 
(a ■ Vu*", u™") = 0, we have 

F(0-e= ^L(/^(u™) + D(a)):Z}(u-) 
(4.5) + 4 |D(u'") + D(a)|f-2(D(u'") + Z}(a)) : D(u'") 

+ (u"^ ■ Va, u"*) + (a ■ Va, u*^) . 

By lemmas 12.11 and 13.21 and Holder and Young inequalities, we obtain 
the following estimates (for some small positive numbers £j and some 
constants C^., c, which may depend on m): 



, u*" ■ Va, u™) I = I (u™ ■ Vu™, a) 
(4.7) |(a-Va,u"^)| <e2|u"^|?,, + C,,; 

|p-2 

>c,/jD(u-)r >ciiu 

then 



(|/^(u-) + D(a)r' (Z}(u™) + D(a)) - |D(a)|^'-2D(a)) : D(u^ 

1)P' 





> 




|u™| 


p + 

l,p 


(4.8) 


> 


Cl 




IP 

l,p 




> 


Cl 


lu^l 


IP 

l,p 




> 


Cl 




IP 

l,p 



|D(u-) + Z^(a)r' (D(u™) + Z}(a)) : D(u^ 
/jD(a)r^D(a):D(u™), 



|^(^)|P)(P-)/. (Jjz^(u'")|.)^/- 

-mlP I . 
^ ll,p ^£3 1 



(4.9) 

(2;(u-) + D{a)) : D(u-) = \D{u^)\^ + D{a) : D(u-) 
> i |D(u-)p - i |Z)(a)p > ci|u-|?,2 - I L |/^(a)p. 
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Then, taking Ei, i = 1,2, 3, suficiently small, and noticing that ju"*!^ ^ > 
c|u™|^ 2 ^ Cil^p (notice that |^| = |u'"|i^2 is a norm in M™), from ( 14 .Sp - 
f lt:9|) we get 

F(O-e>Ci|e|'-C2>0 
for all G M*" such that |,^| > a/ C2/C1, for some positive constants 

Cl, C2. 

Next, we notice that |u'^|ip is uniformly bounded with respect to 
m. Indeed, multiplying fl4.4p by C,j and summing in j from 1 to m, we 
obtain, as in fl4.5p . 

l/jD(u-)p + i^D(u™):D(a) 

+ (u™ ■ Va, u™) + (a ■ Va, u"*) = 0, 

and then, proceeding with similar estimates to obtain (14.61) . (14.81) and 
(14. 9p . we arrive alQ 

(4.10) \D{un\' + \nni, < c, 

for some constant c. Thus, there exists a subsequence of {u"*}, which 
we still shall denote by {u™}, and a vector field u G I'o''''(n) such that* 

inPo^'^(fiT) 
^ ^ u'^^u inL«(fiT) 

when m — )■ 00, where g > 1 is any number less than the critical Sobolev 
exponent p* := = j^, if ri = 3 and p < 3, and 1 < g < 00 is 
arbitrary, if p > (n = 2, 3). In particular, 1 < g < 00 is arbitrary for 
n = 2, since p > 2. 

Now we want to pass to the limit in (14. 4 p when m — t- 00 and obtain 
it with u in place of u"* and with any ip G in place of ip^ . We 

begin by defining the operators 

A{w) = -div{|D(w) + D(a)|P-2(D(w) + D(a))}, 

and 

C(w) = -ldiv{D(w) + D(a)}, 

for w G Vl'^{Q). More precisely, A and C are operators from Vl'^{Q) 
into Vl'^{Qy, defined by 

(A(w),(^)= / \Diw) + D{a)r\D{w) + D{a)) : Dicp) 

and 

(C(w),(^)= / (D(w) + D(a)):D(c^). 
*Here we write explicitly VIt, instead of n, for future reference. 
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Notice that D{w) + D{sl) G Lp' (il) because p > 2 ^ p' < p and = fi^ 
is a bounded domain. We also write 

B(w) = — (w ■ Vw + w ■ Va + a ■ Vw + a ■ Va) . 

So we want to show that {A{u) + C(u), cp) = (-B(u), ip) for every <p in 
©(fi), or, equivalently, for every if in Vq^{Q). 
By (^i^, we have 

(4.12) (A(u-) + C(u™), cp) = (5(u-), cp) 

for all (p e Vl'P{fl) and all m = 1, 2, ■ • • . 

Since |u'^|i^p is uniformly bounded, by Holder inequality we have that 
{^(u'")} is a bounded sequence in Vl'^{Qy, so there is a x ^ ^o'^(^)' 
and a further subsequence {u"*} such that 

(4.13) {A{un,cp) {x,<p) 
for all ip G Next, we show that 

(4.14) (i?(u'"), cp) {B{u), cp), ^pe V'/iSl). 

Fisrt we notice that 

Ku*"- Vu'",£^) - (u- Vu, ip)\ 
= [((u^-u)- Vu"^,cp) + (u-V(u™-u),c^)| 
= |((u™-u) ■ Vu™,cp) - (u- V(^,u'"-u)| 
< llu- - ullJI Vu-||p||c^||, + ||u||,|| V(^||p||u- - u||, 0, 

where q is large enough such that - + - < 1 and less than p* := 

"2 q p — ^ n—p 

if p < n. Notice that if p < n then = 3 (n = 2, 3 in this paper and 
p > 2) and, since p > 2, we have ^ + ^ < 5/6. Similarly, and more 
easily, we also have 

|(u™-Va,cp)-(u-Va,(^)|^0 

and 

|(a-Vu™,(^)-(a-Vu,c^)| -^0. 

Thus we have shown (14.141) . From (14. lip and the fact that p > 2 and 
Q = Qt is bounded, we also have lim(C(u''"), = {C{u.),ip) for all 
ip e Vl'^i^). Then, from flCTjl - flTO]) . we have x + C{u) = B{u) 
in Vl'^{^iy. Then, to conclude the proof, it remains to show that x = 
A{u). To see this, it is enough now to show that {A{u^), u™) converges 
to (x, u), since, by Lemma ^72\ the operator A is monotone. Indeed, 
we have the following classical argument for monotone operators. From 
(A(u'") - A{w), u"^ - w) > 0, i.e. 

(A(u™),u'") - (A(u"^), w) - (A(w),u™) + (A(w), w) > 0, 

if (y4(u™'), u™) converges to (x? u) then, by (14. lip and (14.131) . we can 
take the limit in this inequality when m — )■ oo and obtain {x — Aw, u — 
w) > 0, for all w G Vq ^{^1). Now replacing w by u — Aw, for A G 
M^, we arrive at {x — ^(u — Aw),w) > for all w G V^'^^Q) and 
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all A G M^. Then the desired result follows, once one shows that 
limA->o+(^(u — Aw),w) = (y4(u),w). Here, we can show this using 
the Lebesgue's dominated convergence theorem, since the integrand in 
(^(u — Aw),w) is dominated, for any A G (0,1), by some constant 
times the function {\D{u)\p-^ + \D{w)\p-'^ + \D{a)\P-^)\D{w)\, which 
belongs to L^{Q). 

To show that (A(u™),u") = (5(u™),u") - (C(u"),u'") converges 
to (x, u) which is equal to {B{u), u) — (C(u), u), we write 

(S(u"),u™) - (5(u),u) 
= [(S(u™), u™) - (5(u), u™)] + (S(u), u'" - u) 

and notice that the two last terms converge to zero, when m — )■ oo, by 
the estimates above we used to obtain fl4.14p . It is easy to see, using 
again (14. lip and the fact that p > 2 and Q = Qt is bounded, that we 
have also lim(C(u'"), u™) = (C(u), u). □ 

Next, given any t > 0, we show that the solution of (14. ip is 
uniformly bounded in I)l'^{Qt), with respect to T, for T > t + 1. Pro- 
ceeding similarly to we introduce the function 

(4-15) y{t) = i |u^|^_,^^^ + |u^|^^^^^^ , t > 0, T > t + 1. 

In the sequel we write U"^ = u and often u + a = v. Multiplyng 
(I4.ip i by u and integrating by parts, using that u\dQ = 0, we have 
(4.16) 

14|D(u)|V4|D(v)r^D(v):D(u) 
= — ^ D{a) : -D(u) + (u ■ Vu ■ a — a ■ Vu ■ u — a ■ Va ■ u) 
+ J,^ (n . (D(v)n) + u ■ (|D(v)p-2D(v)n)) 
- /r, (i|uP(u ■ n) + (u ■ a)(u ■ n) + (u ■ n)P), 

where Tt = S(t) U S(— t). First we estimate the 'interior' integrals 
Jj^^ ■ ■ ■ . Using Young inequality and Lemma 12. 1[ we get 

(4.17) - /^^ D{E) : D{u) < e J^^ \D{u)\' + C„ 

and 
(4.18) 

f^^ (u ■ Vu • a — a ■ Vu ■ u — a ■ Va ■ u) 

< \^h,P,n. ( In, |ar>ry^'' + (a ■ Vu ■ u - a ■ Va • u) 

where e > is fixed below. Besides, proceeding as in (14.81) . we get 
,4ig) /„.|Z)(v)rZ)(v):C(u) 
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and 

(4.20) |4 |D(a)r-^D(a) : D{u)\ < e \n\l^^^^ + C,t. 

Then, from fl4.16p - fl4.19p and taking e ^ 1, we obtain 

(4.21) y{t)<cit + I, 
where 

(A oo\ ^ = k [1^" ■ + " ■ \DMr'Div)n - l\u\\u ■ n) 

^ ^ -(u-a)(u-n)-(u-n)P]. 

Now the idea is to control the boundary integral I by the interior 
integral y{t), but if for instance one tries to apply the trace theorem 
then higher order derivatives arise. To achieve that purpose we use 
the clever idea given in ^15j for the case p = 2, that is, to integrate 
I = I[t) from ?7 — 1 to ?7, for 77 > 1, or better, integrate the estimate 
f l4.2ip . Thus we introduce the function 

(4.23) z{7]) = r y{t)dt. 

Jr]-1 

Notice that since y is a, nondecreasing function we have y{7] — 1) < 
^iv) ^ y{v) for S'll ^7 > 1; thus estimating y is the same as estimating 
z. Another interesting feature of the function z is that 

(4.24) z'{7]) = y{r]) - y{r, - 1) = ^|u|2 ^ ^ + 

Then if we estimate I(t)dt in terms of |u|^ ^ ^ and f;|u|^ 2,o^_i ^, 
in the end, in virtue of fl4.2ip . we shall obtain a estimate for z{r]) in 
terms of z'{7]). Then we shall use Lemma [371] to get the desired estimate 
for z{ri). Let's do the details. 

By flOTD and and the fact that J^_^ J^^ ■ = J^^_^ ^■,we have 

n 1 

(4.25) z{r])= y{t)dt<ci7] + -h + l2 + h + h + h, 

Jr)-1 ^ 



where 



^1 = L„_i„u-^(v)n 

/, = ^^_^;^u.|D(v)r2D(v)n 

^3 = J|u|2(u-n) 

^4 = -/f5^_^^(u-a)(u-n) 

^5 = -k..^^nn-n). 
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Using Holder inequality, Lemma I2.1l ii). Poincare inequality (Lemma 
I3.2l iii)) and Young inequality, we have 

< c \D{u)r'\u\ + \D{^)r'\u\ 

^ I I ip— 1 II II I |rj|P^l II II 

- ^ U"ll,p,n^_l,^ll"llp,!^r,-l,r, + l^ll,p,n„_l,^ll"llp,f^»,-l,r, 

SO 

(4.26) \h\<c{z\^)^2!{^)'l^). 
Analogously, 

(4.27) \h\<c[z!{i])^z\i]fl''). 
Regarding J3 and l^, using Sobolev embedding, we get 

I/3I + I/4I < L,_,J|u|' + /o„.,,c|u|2 
(4-28) < c(|u|?,p,^^_^^ + |u|?,^,^^_^J 

= c {z'{rifl^ + z'ijifl'P) . 

To estimate Is, we use Lemma Let w be a vector field in WQ'^{Qr]-i,rj) 
such that V ■ w = u ■ n and |w|ipQ^_^^ < c|u|p^f7^_j^, where c is some 
constant, independent of w and u. Then, using the equation fl4.1l) i . we 
can write 

h = -L V{u -11)= L VPw 

^ i-^ri — 1.71 '^i'rj — 1.7? 



I^(v)|-^^(v) -DM + J^^^_Ju . Vu) . w 



+ ■ Va) ■ w + ■ Vu) • w + /^^_^ ,^(a ■ Va) ■ w. 

Thus, proceeding with similar estimates to those used to obtain (14.271) - 
(I4.28p . we arrive at 

(4.29) I/5I < c {z'{r]) + z'irjy/^ + z\r]fl^ + z'ir^f''^) . 
From (1125])- (1129]), we have 

(4.30) z{t]) <cir] + m {z\r])) , 

for all > 1, with ^(r) = C2 (r + t^'^ + t^/p + r^/P + r^/p) . Now, from 
(I4.10p and the weak convergence (I4.11l) i . we have y{T) < c for some 
constant c (independent of T), so by z{T) < y{T) and by assuming 
that ci > c, without loss of generality, we have 

(4.31) z{T) < ciT. 
Next, let C3 > be a constant such that 

(4.32) 2ci + C3 > 2^(2ci). 

Then, by (I4.30p - (l4.32p . we have the conditions of Lemma l3.1l i) satisfied, 
with ip{ri) = 2ci?7+C3, 5 = 1/2, to = 1 and T > 1 (arbitrary). Therefore, 
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< 2cirj + C3 for all 77 > 1, and hence, since yij] — 1) < zirf)^ we 
conclude that there are (new) constants Ci,C2 such that 

(4.33) y{t) := ^|u^|?,2,n. + l^^l W < cit + C2, 

for alH > and T > t + 1. 

Having the estimate f l4.33p . we complete now the proof of our main 
result. 

Proof of Theorem 12.21 Let u*^ be the solution of (14. ip in fi^, k = 
3, 4, ■ ■ ■ , whose existence is assured by Proposition 14. ![ and set = 
in fl/Qk. By (I4.33p . for each j = 2, 3, ■ ■ ■ , the sequence {u''}k>j+i 
is weakly compact in W^'^^flj), thus, by a diagonalization process we 
obtain a subsequence, which we also denote by {u'^}, and an u in 
W^^^m such that 

(4.34) 

for any t > 0, where g > 1 is arbitrary, if p > n, and less than 
p* := J^, if n = 3 and p < 3. (Cf. (jUID). Besides, by the 

estimate (I4.33P and the fact that u'^ G VQ^iyL), we have that the limit 
u satisfies (I2.2p 9-f l2l2|) t^. Then, to conclude the proof of Theorem 12.21 it 
remains to prove that u satisfies the equation (I2.2p i . in the weak sense 
(12. 4p . Again, we shall use the Browder-Minty method, due to the shear 
dependent viscosity. The idea here is to mimic the proof of Proposition 
14. 1^ paying attention that now VL is not a bounded domain and -D(u) 
is only locally integrable in Vt. This lead us to localize the arguments 
and operators used in that proof, as follows. 

Given 93 G ViVL), letting /cq ^ N such that supp Lp C ^ko-i, we have 



(4.35) / §fc(u'=) : D{^) = / 5(u'=) • yj, 

for all k > kn, where 



§/c(w) = Q + \D{w) + D{ei)r'^ {D{w) + D(a)) 

and 

B(w) = — (w ■ Vw + w ■ Va + a ■ Vw + a ■ Va) . 

Then, we want to pass to the limit in (14.350 when /c — )■ 00 and obtain 
(12. 4p . Let ( : Q — > M+ be a smooth function such that C = 1 in 
supp (p and ^ = in i7 \ Qko and A(^, A(^^k be the operators defined by 

(A^,fcWi,W2) = L Sfc(wi) : D(w2)C 
(A<;Wi,W2) = L S(wi) : D{w2)C, 

on the space 

Vo = w''P{nk,, an n anto) ■= {w g W'^ink,) ■,w = omann an^o} , 



20 G.J. BIAS AND M.M. SANTOS 

where 

§(w) = |D(w) + D{a)f-\D{w) + D(a)). 
Thus, (14.351) becomes 

(4.36) {A^A^'),<fi) = {B{u'),<p) 
and (12.41) becomes 

(4.37) {Mu),cp) = iBiu),cp). 

We notice, as C is a nonnegative function, that ^ is still a monotone 
operator. Besides, {^^ ^(u'^)} is a bounded sequence in V^*, then, up 

to a subsequence, we have A^_fc(u^) ^ xc ^oi some xc iii KT- in 
(I4.14p . we also have 

(4.38) (5(u'=),cp) ^ {B{u),cp). 

Then, by (I4.36p . we obtain (x^, (p) = (-B(u), (p), so it remains to show 
that Xc — ^c('^)- obtain this, from the monotonicity of A^.^' 
enough to prove that (A^_fc(u'^), u^) converges to {x(,u). Indeed, 

(4.39) (Ac,fcU^ u^) - (Ac,,u^ w) - (^^.^w, u'^) + (A^^w, w) > 0, 

for all w G Vo and, by (14.341) . {A(^±w,u^) and (A^^^w, w) tend, respec- 
tively, to (y4^w, u) and (A(^w,w), when k — >■ oo. Then, once we have 
limfc^oo(^C,fc("^)'"'') = (XC'W> we shall have {X(;-A(;{n-\w),w) > 
for all w G Vo and all A > 0, and by Lebesgue's dominated convergence 
theorem, limA_i>o+(^c(u — Aw), w) = (y4^(u),w), hence Xc = ^c(^)- 
Let us show then that \imk^oo{A(^k{u''),u'') = (xc,u). We compute 
(Xc, u) and limfc^oo(^c,fc("^)' '^'^) using directly the equation (I4.1l) i . 
with T = k. Multiplying this equation by and integrating by parts 
in fifcg, we arrive at 

(A^,fcu^u) 

(4.40) = K ■ - K ■ - ^ K " ■ ■ 

-i {K ^ ■ ■ + K " ■ ■ ^0 ' 

where V'' is the pressure function associated with u'^. From (I4.33p . we 
have 

u-D(u*^)-VC+ /" u-D(a) • VC 1 < ^ ^ 0. 




k 



«0 



and that {S(u'^)} is uniformly bounded in L^'(fifc(,), so there is a Xp' ^ 
(fifco) such that 



(4.42) 



^lim f u.§(u'^)- VC= / u-xp'- VC- 
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Similarly to the proof of (14.141) . we also have 

(4.43) lim f B{u'')-Cu= I S(u) ■ (u. 

Next, we show that 

(4.44) / V^Vi-VC^ j Vu-VC, 



«0 



for some further subsequence of — > oo, where, up to a constant, V 
is the pressure function associated with u. For this, it is enough to 
show that there is a P E (flko) {p' = p/{p — 1)) such that "P^ ^ 
V in LP (fifco); {V''} is uniformly bounded in i^ko)- Let us 
assume, without loss of generality, J^^ V^dx = 0. Writing 

g = \vr'-'v'-\n,,\-' f \vr'-'v'dx, 



by Lemma 13.41 there exist a constant c (independent of k) and a vector 
field V e WQ'^iQko) such that 



(4.45) 



V- V = 



1 

yk 1 1 p-i 



l,p,Qfe(, - ^ II' llp',Ofcg • 



Notice that J^^jdx = 0, g e L^i^k,) and \\g\\p,n,,^ < 2\\P''\\;r^^^. 

Then, 

(4.46) 

= gV'dx + |fi,J-i (4^^ \vf-'v') V'dx 

where, for the last iguality, we used equation (I4.ip i . Using again (I4.45P 
and previous estimates, it follows that 

< 

Therefore, 



(^•^'') n" fcii "ii k\\¥^ 



IIp'^II < c llu'^IL „ + ||a|L o < C, 

as we wished. 

From f O0|) -f lO4|) . we obtain 
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Now, replacing u by u'^ in fl4.4Up . we have 
(4.49) 

-i(4„"-^(^')-VC + L,^^n.Z)(a).VC), 

and taking the limit when — )■ oo in the right hand side here, analo- 
gously to the steps we did to obtain fl4.48p . we get the right hand side 
of flAj) . i.e. 

" ■ ^("') ■ VC + 4^ u . D(a) . VC) } 
= k^ B{n) . (Cu) - /^^^ Vu-VC-l J^^^ u . X,' ■ VC . 

Then, combining fl4.48p and (14.491) . we have limjt^oo(^c.fc'^'^' ~ 
{X(, u), and thus conclude the proof of Theorem 12.21 □ 

Next, we make some remarks and prove two additional results, one 
on the rate of dissipation of energy of the solution obtained for problem 
(12. ip and another on the uniqueness of solution. 

Remark 2. Dropping the convective term v ■ Vv in (I2.ip i . we obtain 
the Ladyzhenskaya-Solonnikov problem for Stokes' system with a power 
law. The solution of this problem can be obtained as in the proof of 
Theorem \2.^ with obviously much less computations. 

The solution of problem (12. ip has energy dissipation uniformly dis- 
tributed along the outlets. More precisely, we have the following result, 
which generalizes Theorem 3.2 in [15j for power law shear thickening 
fluids. 

Proposition 4.2. Let \ be a solution of problem (12. ip . with p > 2, 
obtained by the proof of Theorem \2.B . Then there exists a constant k, 
such that 

(4.50) [ |Vv|P <K, Vt > 1, 

Jni,t~i,t 

where i = 1,2. 

Proof. Let u = v — a. By the proof of Theorem 12. 2[ u is the weak limit 
in Wi^^{Q) of a sequence {u'^}^^, where u'^ is a solution of problem 
(14. ip with T = k. Now, for r > max{t, 2} we define the function 

rv 

Zriv) = / yT{t)dt, r]>l, 

Jr)-l 

where 

yJt) := - \u^\, + u 1 o 
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(see Section [2] for the definition of fli,r-t,T+t) ■ Similarly to the proof of 
( I4.33p . it is possible to show that 

Zriv) < Viji) , Vt] G [l,r] , 
where ip{ri) = C2?7 + C3, for some constants 02,03. Since 

/•3/2 1-3/2 

yril/2)= / yril/2)dt< / y^it) dt = Zr{3/2) < ip{3/2) = c , 

Jl/2 Jl/2 

we have 

rT+1/2 



/ / \Vm'\' <yr{l/2)<c 

Jt-I/2 J^i 



-1/2 

and, consequently, by the weak convergence of to u, we also have 

f.T+1/2 

|Vur<c. 



T-1/2 JS, 

which is fl4.50p with u in place of v. Since, by Lemma I2.H the vector 
field a also satisfies this property, this end the proof of Proposition 14.21 

□ 

In p. 1437], Marusic-Paloka observes the difficult of obtaining 
uniqueness results for Navier-Stokes system with a power law. In par- 
ticular, this is an open question even in bounded domains. We can 
prove an uniqueness result for problem (12. ip under some conditions. 



which we specify precisely in Theorem 14.41 below. One of these con- 
ditons is motivated by Proposition 14.21 and another, by the following 
propostion, which was inspired by the solution of Leray problem given 
by Marusic-Paloka; cf. [HI Lemma 4.2/(4.24)]. 

Proposition 4.3. For i either equal to 1 or 2, let v = {vi, ■ ■ ■ ,Vn) be 
a divergence free vector field in Wl^^{Qi) vanishing on dVti and having 
property (I4.50p . If for some j G {1, . . . , n} and some positive number c, 



(4.51) 



for all ) eVLi, then there is a constant C such that 

(4.52) (w Vv,w)n.,, < Ck\\ |Z}(v)|(^~2)/2Z}(w)|||^^^^, 
for all w G Vjf^{ni) andt>0. 

Proof. Denote Q = Qi. By Holder inequality and (I4.50p . we obtain 

(4.53) |(w- Vv,w)n,_,,J < C4w\\l,^^^_^^^. 

By Sobolev embedding and Poincare inequality (Lemma 13. 2p . we have 
(4-54) ||w||2p',ni_i,t < <^|w|i,,.,n,_i,, , 
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for any r G (1,2) such that 2p' < Now, by Korn inequahty 

(Lemma 13. 2p . Holder inequahty and (I4.5ip . we obtain 

\'^\irnt-it~ Iq iVwl*" 

< ciiiz^(v)|(^-2)/2D(w)r,,^^_^,^(/„'^ s-'-Ti^ds)('-^y'. 

Then, chosen r < it fohows that 

(4.55) |w|i,,,^,.,, < C|||D(v)|(^-2)/2D(w)|h,n._,,. 
Thus, from (14.541) and (I4.55p . we have 

l|w|bp',f.._., <C|||D(v)|(^-2)/2D(w)|h,n._., 
and from (I4.53p . we get 

(4.56) |(w Vv,w)^,_,J <CMI |/^(v)|(^-2)/2Z)(w)|||^^_^^^. 

FinaUy, writing fi^ as a finite union of domains Qt-j^i^t-j, 
j = 0, ■ ■ ■ ,m < oo, and adding inequahty (I4.56P with Qt-j-i,t-j in 
place of t with respect to j, we obtain (I4.52p . □ 

Remark 3. An example of a solution satisfying property (14.51 1) when 
Qi is a straight outlet (i.e. the cross sections S(a;i) are constant, with 
respect to xi) is the Poiseuille flow in fij. See [TTl §3]. 

We now state and prove our uniqueness result. 

Theorem 4.4. Let k > be sufficiently small and I be some positive 
number. Then there is no more than one weak solution of problem (12. ip 
inWf^^{Q,) and satisfying fl4.50p and property (I4.52p m fi, i.e. for some 
constant C , 

(4.57) (w- Vv,w)^, < Ck\\ |D(v)|(^-2)/2d(w)|||^^, 
for all w G V]f^{n) andt>0. 

Proof. Let vi and V2 be solutions of (12. ip satisfying the assumptions 
in Theorem 14. 4[ Denote w = Vi — v. Then 

-div{Z}(w) + \D{w) + D{^r2)\f-^[D{w) + /^(va)] - \D{^r|)\P-^D{Y2)} 
+W • VW + W • VV2 + V2 • Vw + V (Pi - V2) = 0, 

where Pi, 7^2 G W^ioc(^)- Multiplying this equation by w and integrat- 
ing by parts over Qt, similarly to derivation of (I4.16p . we obtain 

4 {\Diw) + D{^r2)r^[Diw) + D(V2)] - \Div2)r^Div2)} : Diw) 

= -(W ■ VV2, W)f7, - / , 

where 

^ = -/an.¥(w-n + V2-n-(Pi-P2)(w-n) 

+ /,^^ w ■ {\Diw) + D(v2)|^-2[D(w) + D(V2)] - |D(v2)|^-2Z^(v2)}n . 
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Estimating some terms in the above equation by using Lema 13.21 and 
assumption f l4.57p . it follows that 



Cl / iVwI^'+Cs / |/^(v2)r'|D(w)^ < C3|«:| / \Di^^2)r'\D{w)\'+I. 
Jfit Jfit Jfit 

for some positive constants Ci, 02,03. Thus, if \k\ < Ci/c2, we have 

yit) ■■= Hl^,nMt) < ^ 

(for some constant c). Next, integrating y{t) from r^ — 1 to r/, r/ > 1, and 
proceeding similarly to the proof of fl4.3Up . but using f l4.50p instead of 
Lemma [2.1l ii). we obtain 

Z{7]) < {z'{7])) , 

with \E'(r) = r + T"^/^ + r^/*' + r^^^. Now suppose z is not identically 
zero. Then, by Lema 13. we have 

lim z{t) = 00. 

Besides, since for r > ri (for some ri > 0) 

by Lemma 13.11 again, we also have 

liminft_>.oo e~*2;(t) > 0, if j9 > 3 

liminft^oot"^/^^"^^^(t) > 0, if 

This contradicts (12. ip ^. Therefore, z = and so, vi = V2. □ 

Remark 4. By tracking all the estimates we did to obtain fl4.33p . simi- 
larly, to obtain fl4.5Up . we can see that the constant k in fl4.50p depends 
on the flux a so that k = 0{\a\'^), for some positive number 7. In 
particular, k, tends to zero when a tends to zero. 

Remark 5. For an example where condition fl4.57p is accomplished, 
see da p.l437/§4.2]. 

Remark 6. Regarding the Stokes system with a power law, i.e. system 
(11. 2p discarding (v ■ V)v, we have uniqueness of solution for the corre- 
sponding Ladyzhenskaya-Solonnikov problem for any flux a, as occurs 
in the case p = 2 [TF, Corollary 2.1, p. 739]. 
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